Perturbative Analysis of Dynamical Localisation 
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Abstract. In this paper we extend previous results on convergent pertur- 
bative solutions of the Schrodinger equation of a class of periodically time- 
dependent two-level systems. The situation treated here is particularly suited 
for the investigation of two-level systems exhibiting the phenomenon of (ap- 
proximate) dynamical localisation. We also present a convergent perturbative 
expansion for the secular frequency and discuss in detail the particular case 
of monochromatic interactions (ac-dc fields), providing a complete pertur- 
bative solution for that case. Our method is based on a "renormalisation" 
procedure, which we develop in a more systematic way here. For being free 
of secular terms and uniformly convergent in time, our expansions allow 
a rigorous study of the long-time behaviour of such systems and are also 
well-suited for numerical computations, as we briefly discuss, leading to very 
accurate calculations of quantities like transition probabilities for very long 
times compared to the cycles of the external field. 



1 General Description and Previous Results 

The study of periodically or quasi-periodically time-dependent two-level systems 
is of basic importance for many physical applications, ranging from condensed 
matter physics to quantum optics, as in problems of the theory of spin resonance, 
in problems of quantum tunnelling or in the semi-classical theory of the laser. 
They can be used, for instance, to describe the behaviour of a spin 1/2 system 
in a time-dependent magnetic field, in which case the corresponding Schrodinger 
equation takes the form (we adopt h = 1) 

i3t# = H(t)V, with H(t) = --B(t)-a, (1.1) 

where *(t) = e C 2 , B(t) = (B^t), B 2 (t), B 3 (t)) and a = fa, a 2) <7 3 ) are 

the Pauli matrices. 

Systems like this have been analysed by many authors in various approxi- 
mations, as in the pioneering works of Rabi § , of Bloch and Siegert [0] and of 
Autler and Townes (To) (see also || ^| for more recent dis cussions) . We should 



remark, however, that the interest in the solutions of (1.1) is not restricted to 



the investigation of quantum systems. As first pointed by Feynman, Vern on a nd 



Hellwarth pll (see also the recent discussion in |4|), the quantum system (1.1) is 
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equivalent to the classical Hamiltonian system describing a classical gyromagnet 
precessing in a magnetic field: -^S = —B{t) x S, where S is a unit vector. 

Of particular interest is the situation where the Schrodinger equation takes 
the form 

»at*(t) = fli(t)*(i), with JT x (t) := ea 3 - f(t)a lt (1.2) 

where /(<) is a function of time t and e € K is constant. By a time-independent 
unitary transformation, representing a rotation of tt/2 around the 2-axis, we get 
the equivalent system 

idMt) = H 2 (t)$(t), with H 2 (t) := eat + f(t)a 3 , (1.3) 

where <&(*) := exp (-i7tct 2 /4) $(t) and H 2 (t) := exp (-z7rcr 2 /4) -ffi(i) exp (i7rcr 2 /4). 



One can either interpret the system (1.2) as describing a spin 1/2 system 



as (1.1) under a magnetic field B — (2f(t), 0, — 2e), or as a system with an un- 
perturbed diagonal Hamiltonian Hq := ea 3l representing a two-level system with 
energy levels ±e, subjected to a time-dependent perturbation Hj(t) := —f(t)ai, 
inducing a time-depending transition between the unperturbed eigenstates of Hq ■ 



The equivalent system (1.3), in turn, represents either a spin 1/2 system as < \l . l| ) 



under a magnetic field B = (— 2e, 0, — 2/(t)), or a two-level system composed by 
two uncoupled (for e = 0) orthogonal time-dependent states exp ^— i J* /(t)g?t^ (J) 

and exp (+i J Q ' /(t)gJt^ (°), subjected to a constant perturbation eci inducing a 
transition between them. 

To explain the purpose of the present paper, we have to describe some of 
our previous results. In 0] and (2| we studied the system described by ( |l.2j ) or 
( pT^ ) in the situation where / is a periodic or quasi-periodic function of time and 
e is "small". It is well know that the usual perturbative approach, based, f.L, on 
the Dyson series, leads to difficulties involving secular terms (i.e., polynomials 
in t that appear order by order in perturbation theory and spoil the uniform 
convergence (in t) of the perturbative series) and, for quasi-periodic interactions, 
small denominators. This last problem is typical of perturbative approximations for 
solutions of differential equations with quasi-periodic coefficients and is well-known 
as one of the main sources of problems in the mathematically precise treatment of 
such equations. 

In [jjj and a special perturbative expansion (power series expansion in e) 
was developed, whose main virtue is to be free of secular terms. The algorithm 
employed involves an inductive "renormalization" of a sort of effective field intro- 
duced through an exponential Ansatz (the function g, to be introduced below). For 
the sake of the reader we will shortly recall our method of elimination of secular 
terms in Section ||. In the general case where / is quasi-periodic, it was established 
in that the coefficients of the expansion are also well-defined quasi-periodic 
functions of time but, due mainly to the presence of small denominators, we were 
not able to prove convergence of our e-expansion. Actually, a convergent power 
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expansion in e is not expected without further assumptions (for a detailed analysis 
of these issues in related systems, see Jl2[). 

Less problematic is the situation where / is a periodic function, when the 
obstacle represented by the small denominators is naturally absent. In we 
showed how the difficulties analysed in jl| can be circumvented in the case of 
periodic / and we were able to establish the convergence of our perturbative e- 
expansion uniformly in t G R. 

As discussed in ||], our method not only recovers th e F lo quet form of the 



solution of the time-depending Schrodinger equation (see ( |l.8| )-( |l.9| ) below), but 
also allows the computation of the secular frequency and of the Fourier coefficients 
in terms of explicit convergent e-expansions, what constitutes a feature of our 
algorithm, compared to other expansion methods. 

Due to the technical difficulties involved, we restricted our analysis in 
to two classes of periodic functions, namely those satisfying the conditions (I) 
or (II) presented below (see [|]). Our purpose in the present paper is to extend 
the results of [|j to an additional class of periodic functions. The inclusion of 
this additional class leads to a essentially complete perturbative solution for some 
simple periodic functions, as /(£) = F + <pcos(ujt), representing the important 
case of a monochromatic interaction (also known as ac-dc fields). 

The situation we treat here is also relevant for the rigorous discussion of 
the phenomenon of dynamical localisation, also known (less properly) as coherent 
destruction of tunnelling. This phenomenon, first pointed in |l4| ], indicates the 
possibility to (approximately) freeze the initial state of a quantum system through 
the action of a suitable external time-dependent interaction. This effect has been 
the object of various recent investigations. In ffl], for instance, a rigorous general 
criterion for the occurrence of dynamical localisation was established and applied 
to interesting situations, like the ac-dc field and the bichromatic field. Some of the 
conclusions of |?J on the ac-dc field are indirectly reproduced in Section ||, below. 
We refer the reader to [Hand [Q for more references on this subject. 

The main result of |g] can be captured in the next theorem, for whose state- 
ment we need a definition we will repeatedly use in this work: for an almost periodic 
function h we denote by M(h) the "mean value" of h, defined as 



/' 

T^oo 2T 



M(h) := lim ^ / h(t)dt. (1.4) 



We remark that the limit in ( |l.4| ) is always well defined for any quasi-periodic func- 
tion h. The mean value M(h) equals the constant term in the Fourier expansion 
of h. Details can be found in || |§. 

Theorem 1.1 Let f be a real T u -periodic function of time (T u :— 2tt/lu with 
lo > 0) whose Fourier decomposition f(t) — X^nez F n e muJt , contains only a finite 
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iber of terms, i.e., the set of integers {n £ Z\ F n ^ 0} is a finite set. Let 

m = ( i*® ) = mm = u{t ' ° )m (L5) 



be the solution of the Schrodinger equation (l.l ). Consider the two following dis- 
tinct conditions on f : 

(I) M(Q ) ? 0. 

(II) M(Qo) = but M(Qi) ^ 0, where 

q(t):=exp(i^ /(r)dr^ , &>(*) := g(i) 2 = ex p( 2 *^ (1-6) 



QoW / (Qo(r)- 1 - M (Q- 1 )) dr. (1.7) 
Jo 



Then, for each f as above, satisfying condition (I) or (II), there exists a constant 
K > (depending on the Fourier coefficients {F n , n£Z,n/0} and on u) so 
that, for each e with |e| < K , there are £ R and T^ -periodic functions and 
such that the propagator U(t) of ( |-?-4 J can be written as 



U(t) = 



U xl (t) U 12 (t) \ ( U u (t) U 12 (t) 



Un{t) U 22 (t) j \ -U 12 (t) Unit) 



with 

Unit) = e- mt u n it) + e mt u+^t), U l2 {t) = e~ iat u^ft) + e mt u+ (i). (1.8) 

The functions uj^ and uj^ /iaue absolutely and uniformly converging Fourier ex- 
pansions 

u±(t) - ^^(nje^, u ±(t) = X;Wu(») e<wt - (1-9) 

nSZ n<EZ 

Moreover, under the same assumptions, Q and the Fourier coefficients (n) and 
lA^ 2 i n ) can be expressed in terms of absolutely converging power series on e. □ 

Let us now discuss the conditions (I) and (II) of Theorem |1.1| . Writing the 
Fourier decomposition of / as /(<) = F + J2n=i cos(nwi) + sin(nwt) 
the set of all possible functions / with a given J and lj can be identified 

with the parameter space R 2J+1 of all real coefficients i*b, f>[ j, 1 < n < J. 
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The (complex) condition M(Qo) = determines a (2 J) or (2J — l)-dimensional 
subset of 5,7,0)) where condition (II) eventually applies. It is also on this subset 
that the more restrictive condition M(Qq) = M(Qi) = should hold, restricting 
the parameter space of / to a (2 J — 1), (2 J — 2) or (2 J — 3)-dimensional subset, 
if it is non-trivial. One should, therefore, expect that successive conditions like (I) 
and (II) would eventually exhaust completely the set $j tU . 

To illustrate all this, let us consider the simplest example, when / represents 
a monochromatic interaction: f(t) — ip\ cos(u>t)+(fi2 sha.(u>t), with (ipi, tpi) e R 2 . A 
simple computation shows that M(Q ) = e 2n f J (^ a ) , where ip := \ftp\~- V </?f, 
Jo is the Bessel function of first kind and order zero and, in this case, 7/ = ipijw. 
Moreover, M(Qi) — for all (<pi, tp 2 ) € R 2 . (See ||] for details). Hence, condition 
(I) is satisfied for all (ipi, <p 2 ) G K 2 , except in the circles defined by ipa — cux a /2, 
a = 1, 2, ... , where x a if the a-th zero of Jo in R + . Condition (II), however, is 
never fulfilled in this case. To achieve a complete solution we have, therefore, to 



extend Theorem 1.1 to include further conditions beyond (I) and (II), holding on 
the circles cpo = ujx a /2. 

The purposes of this paper are to identify the first condition following (I) 
and (II) , which we call condition (III) , to show that the method of elimination of 
secular terms holds in this case as well and, for periodic interactions, to show that 



the expansion (2.5) converges for |e| sufficiently small, uniformly for t 6 R. As we 
will discuss, this leads to a complete perturbative solution for the monochromatic 
interaction. As we will see, the identification of condition (III) and the application 
of the method of elimination of secular terms to it are highly non-trivial tasks. 



This paper is organised as follows. In Section g we describe the general strat- 
egy employed to eliminate the secular terms and present our main theorems. In 
Section [| we prove some interesting and useful mathematical results on the mean 
value of some quasi-periodic functions, introduce the "renormalisation" operation, 
important to organise our procedures, and introduce some notations we will use 
throughout the paper. In Section |], which is the technically central piece or this 
work, we apply our strategy of elimination of the secular terms to the situation we 
wish to analyse. In Section [5] we apply our results to case of monochromatic inter- 
actions and discuss the issue of dynamical localisation in this case. In Section p| we 
briefly describe some numerical calculations based on our results. Sections |^and 
contain some of the physical applications of our work. In Appendix [X] we treat 
several results used in Section ^, whose proofs unfortunately involve too many 
algebraic computations to be included in the main text. Appendix ^ is dedicated 
to the proof of convergence of our expansions for the periodic case. Appendix ^ 
contains some comments on the Fourier coefficients of the wave functions. Finally, 
Appendix ^ sketches the proof of an identity on Bessel functions we use in Sec- 
tion m 
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1.1 Comments on the Notation 



In this paper, Z + will denote the set of all non-negative integers (zero included) 
and Z* the set of all integers, excluding zero. is the set of all positive integers. 
These notations are also applied to lA and to the real line R. 

Vectors in Z A (or R A ) will be written as v. The operation v ■ u will denote 
the scalar product in 7L A (or R" 4 ), defined as v ■ u :— v\Ui + • • • + vaua- 

For a quasi-periodic function h : R — > C, we write its Fourier decomposition 
as h(t) = J2mei, A Hme v —- ht , where A is some positive integer and u h £ R A (see, 
e.g., jl5| ^6|). The Fourier coefficient Hq will be denoted simply as Hq. 

For m G Z we denote by <m> the following function: 

f ItoI, for m^O , 
<m> := <^ 1 1 ' T . (1.10) 

I 1, for m = 



Beyond the functions Qq and Qi defined in (1.6) and (1.7) we will frequently 
use the following functions 

Q 2 (i) := Qo(t) [ (Qo(r)-M(Q ))dT, (1.11) 
Jo 

Qs(t) := Go(*) / (Qi(r)-M(Q 1 ))dr. (1.12) 
Jo 

Note that, by their definitions, the functions Qi, i = 0, ... ,3, are quasi-periodic 
if / is quasi-periodic. We will have more to say about their properties below. 



2 Elimination of Secular Terms. The Main Re- 
sults 

We recall in this section some of the methods and techniques developed in the 
previous works JjJ and A key result for our method is the theorem below, 
proven in M, which presents the solution of the Schrodinger equation ( [1.3] ) in 
terms of particular solutions of a generalised Riccati equation. 

Theorem 2.1 Let f : R -> R, f e C X (R) and e e R and H g : R -> C, 

5 £ C 1 (R), 6e a particular solution of the generalised Riccati equation 

G - iG 2 - 2ifG + it 2 = 0. (2.1) 

Then, the function <S> : R -> C 2 given &j/$(i) = = C(*)$(0) = C7"(t, 0)$(0), 

w/iere 

i?(t) (1 + ig(0)S(t)) -ieR(t)S(t) 
H/l := I , . I • (2.21 

-ieiZ(t) 5(t) ( 1 - t ,g(0) S(t) 



G 



with 



R(t) := exp (~i J (f( T ) +9( T )) dT ^j and S(t) := R{t)~ 
is a solution of (1.5) with initial value $(0) = f ^ + (oj^ £ C 2 . 



dr, 



□ 



Let us briefly describe some of the ideas leading to Theorem |2.l| and to o the r 
results of jjj. As we saw in [Q, the solutions of the Schrodinger equation Jl.3| ) 
can be studied in terms of the solutions of a particular complex version of Hill's 
equation: 



4>(t) + (if(t) + e 2 + /(i) 2 ) 4>{t) = 0. 



(2.3) 



In fact, a simple computation shows tha t the components satisfy <f>± + (±if + 
e 2 + f 2 )4>+ = 0. If we attempt to solve ( |2.3|) using the Ansatz 



(2.4) 



^(t)=exp (-ij (f(r)+g(T))dT 



it follows that g has to satisfy the generalised Riccati equation (2.1). We then try 
to find solutions for g in terms of a power expansion in e (vanishing for e — 0) like 



g(t) = q(t)J2v n (t)e* 



(2.5) 



71 = 1 



where the function q was defined in (1.6) and i s of central importance in this work. 
The heuristic idea behind the Ansatze (2^4) and (2.5) is the following. For 

e = a solution for (2.3) is given by exp ^— i j* /(r)dr^ . Thus, in (2.4) and (2.5) 

we are searching for solutions in terms of an "effective external field" of the form 
/ + g, with g given in terms of a convergent power series expansion in e, vanishing 
for e = 0. A solution of the form (2.4) leads to one of the two independent solutions 



of (2.3). The full solution of (1.3) in terms of solutions of the generalised Riccati 
equation (2.1) is that described in Theorem 2.1 (see the discussion of [Q). 



We then proceed inserting (2.5) into (2.1 



). The result is a set of recursive 
first order linear differential equations for the functions v n which can be easily 
integrated. The solutions of these equations are 



vi(t) = Kxq(t) 



V2(t) = q(t) 

v n (t) = q(t) 



i f (^iQo(t) - Qo^r 1 ) dr + k 2 
Jo 

i I / v p (T)v n _p(T) dr ) + K n 
. \ Jo p=i J 



(2.6) 
(2.7) 



, for n > 3 , (2. 
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where the k„'s above, n = 1, 2, . . . , are arbitrary integration constants. Defining 
Za(t) := «?CoW- 2o(i) _1 , V*) := X] w ? ®' n - 3 ' 



we can write (2.6)-(2 



as 



= niq(t), v n {t) = iq(t) / T n {r)dr + K n q(i), n>2. 

Jo 

Observe that, in particular, we could just set all the k„'s equal to zero. However, 
this is not a clever choice, since it would result in polynomial terms on t (the 



so-called Secular Terms) for the series expansion (2.5) of g. This, of course, would 
restrict the convergence of the series just for small values of time. As noticed in 
[Q, there is a choice of the constants K n for which one can eliminate completely 
all the polynomial terms on t that would eventually appear in g. This procedure, 
which we call the Elimination of Secular Terms, will be briefly described now. 

First of all, assuming that the function / jl — > R is quasi-periodic, it was 
proven in Appendix B of ffl that q, defined in (2.5), is also quasi-periodic. Hence, 



vi in (2.6) is quasi-periodic. The same is true for the integrand Ti which appears 
in i>2, equation ( |2.7| ). Recalling that X2 depends on the free integration constant 
Ki, the key idea is to fix k\ in such a way that the mean value of Z2 is equal to 
zero, that is M(X 2 ) = M{k\Q — Qg 1 ) = 0. Since Q is a quasi-periodic function, 
it readily follows from this that 



^-"1^ (29) 
_ M(Q ) ' (2J) 

With this choice of K\ one guarantees the absence of a constant term in the Fourier 
expansion of Ti. Since Ti is being integrated in time, this would imply the absence 
of a linear term on t in the final expression for v±. An important remark is that 



(2.£) will only make sense if we assume M(Qq) ^= 0. 

Under this assumption we can now proceed and fix recursively all inte- 
gration constants K m 's by imposing a zero mean value for the integrands X„'s, 
n = 3, 4, . . . , which appear in ( |2.8D . This procedure removes, order by order in e, 
the presence of the secular terms in the series expansion ( |2.5|) for g and recursively 
implies that all functions v n are quasi-periodic. Once all secular terms have been 
removed, one can write the Fourier expansion for the functions v n as 

v n (t) - ]T V^e^\ (2.10) 

provided the sum converges absolutely. It was shown in jj] that this is indeed true. 
The proof of this fact was performed in the following way: first it was shown that 
the Fourier coefficients Qm of the function q satisfies the bound \Qm\ < Qe _x '— ', 
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for some \ > 0- Then, the method of elimination of secular terms described above 
was applied to fix the integration constants n n leading to inductive bounds of 
the form < /Cne^* -5 '*'''— ' for the Fourier coefficients of the function v n , 

where < 6 n < x, for all n = 1,2,.... This exponential decay is enough to proof 
the convergence of the sum in ( [2.1C| ) and to establish by induction the quasi- 
periodicity of all the functions v n . Unfortunately, due to the bad behaviour in n 
of the constants K, n , it was not possible prove the convergence of the e-expansion 



(f2.5|), hence (2.5) has to be seen as a formal quasi-periodic power series solution 



of the generalised Riccati equation (2.1) 



The reason for the bad behaviour of JC n is related to the presence of con- 
volutions and to the small denominators appearing in the recursive relations for 
the coefficients Vnfi . A general discussion of these problems is found in [0 . How- 
ever, in the situation where / is a periodic function, stronger results are possible. 
In 0, where this situation was studied, it was possible to prove the convergence 
of the power series ( [2.5| ) and uniform convergence of the Fourier series involved 
in the computation of the wave functions. Moreover, absolute convergence of the 
e-expansions leading to the secular frequency and to the coefficients of the Fourier 
expansion of the wave functions was also proven. 

All the work done in jy and H was restricted to one of the mutually exclusive 



conditions (I) and (II) of Theorem 1.1. These conditions are consequences of the 
method of elimination of secular terms. Clearly, (I) is vital for ( |2.9| ). When (I) is 
not satisfied we have to apply condition (II). Both cases were studied in [Q and [||, 
where the method of elimination of secular terms has been applied and equivalent 
results concerning the solution g were obtained. In the present work, we apply 
the method of elimination of secular terms to study a more restrictive condition 
than those represented by (I) and (II). Namely, we are concerned here with the 
situation where M(Qo) = and also M(Qi) = 0. 

In this case, the complexity of the calculations involved to find the right 
choice of constants k„'s grows enormously, in contrast with those needed in the 
cases (I) and (II), already studied. The reason for that is quite simple: due to the 
hypothesis M(Qo) = M(Qi) = on e needs to work explicitly with higher order 
terms involved in the expansion (2.5). This will become more clear in Section |J. 



In Section ^ we will discuss an important example where conditions (I) and 
(II) are not satisfied. To solve it, we have to apply the solution (free of secu- 
lar terms) obtained here. We present in Section ^| numerical calculations on this 
particular example and obtained some interesting results. 

We are ready now to state two main theorems of this work. 

Theorem 2.2 Let f : M — > R be a real quasi-periodic function satisfying 
(III) Af(Qo) = M(Qi) = 0, but M(Q 3 ) ^ 0. 



Then, there are constants K n , n > 1, such that all functions v n given in (2.t) 



(2.8) are quasi-periodic. The explicit recursive expressions for the constants n n 
are found in ^.3^-^.3^). □ 
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The proof of this theorem is the main content of Section It states that 
the procedure of elimination of secular terms outlined above also works under 
condition (III). When / is quasi-periodic this does not imply, however, that the 
formal solution Q2.5| ) of the generalised Riccati equation (2.1) converges, since we 
have the same difficulties discussed in detail in JD. 

For periodic /, the situation is different and stronger results can be proven. 

3 imujt k e a rea j periodic function with frequency u>. If Fq = 



Let /(i)=£ m£Z F, 



M(f) — 0, q and Qo = q 2 are also periodic and their spectra of frequencies are 
subsets of {nuj, n € Z}. Following the notation employed in [0, we write the 
Fourier expansions of q and Qo as 



Q (t) = q(t) 2 = 



Q 



(2) gimujt 



(2.11) 



By relatio ns (p^)-(|2.8| ) and with the choice of constants n n mentioned in Theorem 
[2~2| (see ( 4.36D - fl4.38 )), the functions v n are also periodic and their spectra of 
frequencies are also subsets of {nu>, n £ Z}. We write their Fourier expansions as 



(2.12) 



In Appendix ^ we prove the theorem below, which justifies our whole proce- 
dure for the case of periodic interactions and establishes convergence of (|2.5|), 

Theorem 2.3 Let f(t) be as above with Fo — M(f) = and such that condition 
(III) of Theorem \2.^ is satisfied. Moreover, assume that the coefficients Q m and 

(2) 

Qm above satisfy the following: for any x > there is a positive constant Q = 
Q(x) such that 



\Qm\ < Q 



<m> 2 



and 



)(2) 



< Q- 



-x\m\ 

'<m> 2 



(2.13) 



for all m G Z, where the symbol <m> was defined in ( Then, with the 



constants n n fixed as in Theorem 
given in j\2.(\)-(2.8) satisfy 



the Fourier coefficients of the functions v r , 



V, 



(n) 



< 



M (Mi) r 



e -x|m| 

<m> 2 



for all m £ Z and all n > 1, for some positive constants Mq and Mi. As a conse- 
quence, the power series expansion l \2.?\ ), representing a solution of the generalised 
Riccati equation (2.1), converges uniformly for t € R, provided |e| < 1/M±. □ 



Remark 2.4 In ffy, ) was established for f is periodic and represented by a 
finite Fourier series. The condition Fq — M{f) = above is not crucial and can 
be eliminated following the procedure described in Jl/. 
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It follows from this theorem that the main consequences of Theorem 1.1 



are 



valid under condition (III) as well. In particular, the Floquet form (l.S) holds and 
the secular frequency f2 and the Fourier coefficients of Ql-S| ) are analytic functions 
of e for |e| small enough. See Appendix |^ for a comment on this. 



2.1 The Secular Frequency 

A feature of our method is that it allows to present the complete e-expansion for 
the secular frequency (also known as Rabi frequency) associated to the solutions 
of Q-Q (see ([□])). One has (see § § §) 

OO 

n = M(/) + M(g) = F + Y, £ n M (qv n ) . 



By Theorem 2.3 above, this expansion is convergent for |e| small enough. The 
knowledge of the complete expansion is particularly important for the qualita- 
tive investigation of the large-time behaviour of that solutions. After some simple 
calculations using (2.6)-(|2.8D one gets 



= F + ertiM(Qo) + e 2 [i«fM (Q 2 ) - iM{Q x ) + K 2 M(Q )] 

+ e 3 [2 Kl M(Q 3 ) + K 3 M(Q )] + 0(e 4 ) . (2.14) 



As we will show in Corollary 3.2 below, M(Qo) = implies M(Q2) = 0. Hence, 
for case (II), 



il = F -ie 2 M(Q 1 ) + 2K 1 e 3 M(Q 3 )+O(e 4 ), 



and for case (III), 



Q = F + 2n 1 e 3 M(Q 3 ) + 0(e 4 ) 



(2.15) 



(2.16) 



Actually, after fixing k% in Section 4.1, we will see that fl = Fq + 2e 3 \M(Q 3 )\ + 

0(6 4 ). 

In this case (III), if one additionally has Fo = 0, then £1 = 0(e 3 ), a fact first 
pointed in This implies long transition times for certain probability amplitudes, 
a phenomenon known as (approximate) dynamical localisation (see (?]] and other 
references therein). In Sections]^ and [6] we discuss this situation for / describing 
a monochromatic interaction. 



3 Properties of the Mean Value and the Renor- 
malisation Operator 

Let us introduce some notations that will be very useful. Since expressions like 
f{0 Jo d£' 9(£') wm °ft en appear throughout our calculations, we define a short- 
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hand notation 



(/Iff), 



f(0 / dTffCO- 



(3.1) 



Moreover, if (/ | g). is quasi-periodic function on t, then M (/ | g) will denote its 
mean value. We also define 



(f\g\h) ( -.= (f\(g\h) ( ,) = /(o / ^'g(n / #'mo 



(3.2) 



Further compositions like {f\ I/2 | • • • | A)c are defined in analogous way, so that, 
for n > 2, 

(A I/2 I ••• |/n)* := (A I (A I ••• |/„) e ) • 



Note, finally, the trivial fact that 

/o(0(AIAI ••• \fnh = (ZoAIAI ••• \fnh- 



(3.3) 



3.1 Properties of the Mean Value 

The following general results on the mean value of some quasi-periodic functions 
(see definition (1.4)) will be used for many purposes in the present work. 

Proposition 3.1 Let a : R ^ C, 6 : R ^ C and c : R — * C be quasi-periodic func- 
tions with Fourier components denoted by Am, Bm and Cm, Uk & ^ A , respectively. 
We have the following statements: 

1. If M(c) = 0, then the Fourier components of the function h(t) := (b \ c) t are 
given by 



' Cm {^n_ F$ n —m) 



in 1 2 



for all n £ 7L A . Moreover, if M(b) = and C m B- m 
m e then M(h) = 0. ~ ~ 

2. IfM(b) = and M(c) = 0, then M(b\c) = —M (c\b). 

3. IfM(c) = 0, then M (c | c) = 0. 

I IfM(a) = M(c) = M (b\c) = 0, then 



(3.4) 

C-mB, n for all 



M(a\b\c) = -M 



b(t) / a(r) dr 



c(r) dr 



(3.5) 



Moreover, if also M(b\a) = 0, then M (a \ b \ c ) = M (c \ b \ a ) . 
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5. If M(a) = M(b) = M(c) = and M (b | c) = M (c | a) = M (a | b) = 0, 
then 

M (a | 6 | c ) + M (b \ c \ a ) + M (c | a | 6 ) = . (3.6) 

6. If M(a) = 0, it follows from (3) and (5) that M (a | a | a) = 0. □ 



The identity (3.6) is very remarkable. We have noticed that some non-trivial 
relations follow from it, specially if one write ( [3.6j ) in terms of the Fourier coeffi- 
cients of a, b and c. 

Proof. We can demonstrate (1) by explicitly computing the Fourier decomposi- 
tion of h(t). Thus, 



h(t) 



E B n- e 



lit: 



E 



E 



E 



proving the first statement of (1). Now, if M(b) = Bq = and CmB-m = C-mB^ 
for all m € Z^f, then 



E 



m J - J — ni 



Since M(h) — Hq, we completed the proof of (1). To demonstrate (2) we simply 
use the first equality of (3.7) to write 



M (b\c) = Y and M(c\b) = V ^ML 



Changing m — » —m in the second equation above, we get the desired claim. 

Statement (3) is a mere consequence of (2) when we take b = c. Statement (4) 
can be proven using (2). Indeed, let h(t) :— (b \c) t . Since M(h) = M (b | c) = 
and M(a) — 0, by (2), we can write 



M(a\h) = -M(h\a) = -M[h(t) a{r)dr 



= -M 



b(t) / c(t) dr / a(r) dr 



which proves the first claim. To prove the second one, all we need to do is to 
interchange the the roles of c and a in the last equality (note that since M (b \ a) = 
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0, the mean value of (c | b \ a ) t is well defined). Finally, statement (5) can be easily 
obtained as follows. Using fl3.5|), writing b(t) = ^ J"' &(r)dr, using definition ( |l.4| ) 
and integration by parts, one gets 

M(a\b\c) = - r lim ^ (J* b(r) dr^j (J* a(r) dr^j Qf* c(r) dr 



+ M (^(j b(T)dTja{t) c(r)dT 
+M ( ( ( b(r)dA ( f a(T)dr) c(t) 



Now, since M(a) = M(b) = M(c) = 0, the integrals in the first line of ( pTq ) are all 
bounded. Thus, the limit T — > oo is zero because of the division by 2T. Applying 
( |3.5| ) to the remaining terms we obtain (5). ■ 
The following trivial corollary is of crucial importance for some of our calcu- 
lations: 

Corollary 3.2 For M(Q ) = one always has M(Q 2 ) = 0. □ 



Proof. If Af(Qo) = then, by flLllfl , Q 2 (t) = (Q I Qo) t - Hence, from statement 
(3) of Proposition O, it follows that M(Q 2 ) =0. ■ 



3.2 The Renormalisation Operator 

For general quasi-periodic functions a\, . . . , a n , the function (a± | • • • | a n )t, de- 
fined above, is not generally quasi-periodic, since an integration performed on a 
quasi-periodic function with a non-zero mean value would produce a (linear in 
t) secular term, which would eventually become a higher degree polynomial af- 
ter further integrations. We will here describe an operation designed to produce 
a quasi-periodic function out of (ai | • • • | a n )t through interactive subtractions of 
the mean value of the functions being integrated, a procedure we call "renormal- 
isation" due to the analogy to the procedure of perturbative renormalization in 
quantum field theory. We will use this procedure of renormalisation in the following 
sections and here we present its definition and basic properties. 

Let ai, . . . , a n be quasi-periodic functions. We define inductively the renor- 
malisation operator 9\ n acting on (eti | • • • | a n ) by 

9tiai(f) := a x {t), 
$K 2 (ai|a 2 ) t := (o x | 9ti(a 2 ) - Af(3ti(o a ))) t = (<*i I «2 - Af(aa)) t , 
9tn(ai | • • • | a n ) t := (ai | ^„_i(a 2 | • • • | a n ) - Af($H„_i(a 2 | • • • | a„)) ) t , 

for n > 2. We will now prove some elementary facts on 5H„ which will be used 
below. The first important observation is that if a\ , . . . , a n are quasi-periodic 
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functions, then £H„(ai | • • • | a„) is also quasi-periodic. This can be easily seen by 
induction, through the obvious remark that the mean value of 5H„_i(a2 | • • • | a n ) — 
M(*R„_i(a2 | • • • | a„)) is zero. Note also that, trivially 

a a 9\ n (ai\ ■ ■ ■ \a n ) = ^Rn(a ai \ ■ ■ ■ | a„) . (3.8) 

The following proposition is a trivial but useful restatement of the definition 
of the m n 's: 

Proposition 3.3 For all n > 2 the following statement holds: if ax, ■ ■ ■ , a n are 

quasi-periodic functions, then 

9t„(ai | • • • | On) = 9*2 (ai | 5H„-i(a 2 | ■ ■ • | a n )) . 
Consequently, for n > 2, 

Kn(oi | ••• | a n ) = 9* 2 (ai | <K 2 (a 2 I 9*2 ( ■ ■ ■ | ^ 2 K-i | a n ) )••■)) • (3.9) 

□ 

Proof. By the definition of 9t 2 , 

9la(oi|5^i(aa| ••• |a„)) = 

= (ai| £««_i(a2| ••• |o B )-M(5t n _i(a 2 | ••• | a„)) ) 
= 9tn(ffli I • • • I On). m 

Relation ( |3.9| ) shows that the operation 9\ n can be obtained by iteration of 
the operation ![R 2 . One also has the following useful 

Proposition 3.4 For all n > 1 £/ie following statement holds: if at, . . . , a n are 

quasi-periodic functions and a n = 9t 2 (& | c) /or quasi-periodic functions b, c, then 

9tn(oi| ••• |o„_i|a n ) = 9Vi-i(oi| ••• I a n _i | b\c). □ 

Proof. For n = 1, let a x = <K 2 (6 | c). Then 9t x ai = ai = (0 |c), triv- 
ially. For n = 2, let a 2 = 9t 2 (& | c). Then, <H 2 («i | a 2 ) = (a x | a 2 - M(a 2 )) = 
(ai I 9\ 2 (b |c) -M(m 2 {b |c))), but, by definition, 

9ta(ai |6 |c) = (ai I m 2 (b \ c) - M(« 2 (6 | c)) ) 

and the statement holds again. For rt > 2, let a„ = 9^ 2 (6 | c). Then, by induction, 

9Wi(ai I ••• |o„_i |6|c) = 

= (ail 5H„(a 2 | ■•• |a„_i|&|c)-M(9t„(a 2 | ••• |a n _i|6|c)) ) 
= (ai| 9^ n _i(a 2 I • • ■ I a„_i I a„) - M($H„_i(a 2 I ■ ■ • I a„_i I a„) ) 
= %t(ai I • • • I a n _i I a„). ■ 
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If a and bx, . . . , b rn are quasi-periodic, a function like 53fe=i ( a I ^fc ) ma y n °t 
be a sum of quasi-periodic functions, even when (a | Ylk=i &fc) = SkLi ( a I &fc) 
is quasi-periodic, since we are not assuming that M(bk) = for each individual 
fc. This fact notwithstanding, the following simple statement holds and will be 
repeatedly used: 



Proposition 3.5 Let a and bx, . . . , b m be quasi-periodic functions. Then 



<K 2 a 



X> = £*M°I 6 *)- 



fc=i 



fe=i 



Consequently, 



fc=l / fc=l 

/or quasi-periodic functions a\, . . . , a„_i and 61, 



X> = £*»( 



ai 



I On-l I & fe) 



(3.10) 



(3.11) 



□ 



Proof. We have 



fe=i , 



m 

J2(a\h-M(b k )) = ^2 (a I 6. 



£(6 fc -M(6fc)) 



fe=i 



fc=i 



Note, in the third equality, that the mean val ue of — M(bk) is zero and, h ence, 
(a I bk — M(bk)) are quasi-periodic. Relation ( 3.11 ) fol lows from ( 3.E )- ( 3.1C ) . ■ 
The following corollary follows from Propositions 3.4 and 



Corollary 3.6 Let a%, . . . , a„_i, b\, . . . , b m and C\ 

functions. Then, 



X* ax 



a n -i 



^2 ( & * 



. , c m 6e quasi-periodic 

m 

^ 9^+1(01 1 ■ • ■ 1 On-x 1 Pfc 1 Cfc). □ 



fc=l 



fc=l 



Remark 3.7 TTie reader should be warned of the following fact: If ax, ■ ■ ■ , a„, 
61, . . . , 6 m are quasi-periodic functions and (ax \ ■ ■ • | a n ) + (°i I ' ' ' I &m) *s afeo 
quasi-periodic, it is not always true that (ax \ • • ■ | a n ) + (61 | • • • | b m ) equals 

!R»(ai| ••• \ a n )+m m (bx \ ■■■ \ b rn ). 

For a counter-example, take ax(t) = I, a%(t) = e lt , a%(t) = 1 bx(t) = ie lt and 
b 2 (t) = 1. One has (ax | a 2 | a 3 ) + (61 | 6 2 ) = e lt - 1, but fH 3 (ax \ a 2 \ a 3 ) + 
9^2 (61 I £»2 ) = 0. Hence, the renormalisation operations are not additive in this 
sense. 
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3.3 More on the Notation and Some Definitions 



With the shorthand notation introduced in (3.1) and the definition of Q in (1.6), 
we see from (Q, ([Til]) and ( pTl2| ) that 

&(*) - K 2 (Q \Q^) t , Q 2 (t) = %(Qo|Qo) t , fia(*) - %(Qo|Qi) t . 



Below, we will often use the following compact notation 

(i\j) t := (Qi \Qj) t , (i\j) t := (Qi \Q]) t , (i\j\k) t := {Q i \Q j \Q k \, 

etc, for i, j, k = 0, . . . , 3. In other words, we simply use the index n of Q n to 
denote Q„ itself. Moreover, by M (i \j) we will denote the mean value of (i \j) t , 
etc. 

Note that for M(Qq) = one has with this notation Qi(t) = (0 \0) t and 
Q 2 (t) = (0 | 0) t and for M(Qi) = one has 

e 3 (*) = (o |i) t . 

We will often write Q3 this way. 

For M(Q ) — 0, other identities will also be at hand. For instance, one has 

(0 |0 |0) t = (0 1(0 |0) t ,) t = (0 |2) t and (1 1 0) t = (2 | S) , (3.12) 



since M (0 |0) = 0, by item (3) of Proposition 3.1. Relations like these will be 
often employed. 

4 The Case M(Q ) = and M(Q 1 ) = 



In this Section we will prove the Theorem [2.2[ Our interest is to study the situation 
complementary to cases (I) and (II), i.e., the situation where one has condition 

(IIIo) M (fio) = and M{Q X ) = 0. 

To remove the secular terms from g, applying the method described in the previous 
section, we will be forced to add a further restriction to (IIIo), namely the condition 
M (0 1 1) ^ 0. 

Recall that the functions q and Q\ depend primordially on the interaction 
/ (see the definitions given in Q and (fL7|)), so conditions (I| (II) or (III ) 



apply upon the properties of /. As we already saw in Section |lj, the function 
f(t) — tpi cos(wi) + Lp2fiin{uot) only satisfies condition (I) or (IIIo), depending on 
the particular choice of the parameters ip\, ip2- We will have more to say about 
this example latter on in Section ^[ Now, let us work with the expansion for g in 
order to remove all of its secular terms. 

Again, our Ansatz to solve the generalised Riccat i eq u ation ( |2.l|) is ( [2.5| ). 



The explicit solutions for the coefficients v n are given in ( |2.6| ), (2.7) and (2.8). One 
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sees immediately from condition (Mo) that v\ and v-i do not suffer from secular 
terms. Indeed, q in quasi-periodic and since M(Qq~ ) = M(Qo) = 0, we conclude 



that the mean value of the integrand Xi occurring in (2.7) is zero. Therefore, the 



integration occurring in the definition of i>2 in ( |2.7P does not produce a linear term 
in t. These fats imply that v\ and vi are quasi-periodic under (Mo). From these 
considerations we see that the condition M(I n ) = 0, n > 3, becomes recursively 
identical to 2p=i M(v p v n - P ), n > 3, since the v n 's become successively quasi- 
periodic when the interactive procedure is run. 

If t his is ach ieved, i.e., if we succeed in fixing M(X n ) = for all n, we can 
rewrite ( |2.6| )-( |2~8| ) in a "renormalised" form: 

(4.1) 
(4.2) 

n>3, (4.3) 



vi(t) = Kiq(t), 

v 2 (t) = qity 1 (iKjQ 2 (t) - iQi(i) + k 2 Qo(*)) 



v„(t) 



l(t)- l ^J2<R 2 (0| v p v n - p ) t + n n Qo(t) > 



where Qo, Q\ and Q2 were defined in (|l.6|), (1.7) and ( 1.11 ), respectively. 



Let us move on and analyse the third order term. According to (2.8) the 
integrand X3 which appears in the definition of W3 is given by 2viV2- Using (4.1) 
and (4.2) we have 



V1V2 — infQ.2 — iniQi + K1K2QC 



(4.4) 



From Corollary 3.2, we readily see that M(J 3 ) — 2M(viv 2 ) = 0. This means that 
the integrand which appears in U3 does not have a constant term in its Fourier 
expansion. Hence, U3 is quasi-periodic. 

Until now we have verified the absence of secular terms in the series expansion 
of g up to order three in e. As we shall see next, for the same to be true up to 
order four, we have to make a especial choice for the value of the constant K\. 



4.1 The Absence of Secular Terms in t> 4 . Fixing k x 



As one sees from (2.8), the integrand in V4 is X4 :— 2viv$ + v\. Since fi,^,"^ 
are quasi-periodic, the mean value of Z4 is well defined. Let us explicitly evaluate 
li using (|j|)-@: 



X4 



2kiQ (0 i / 2vi{t)v 2 {t) dr + k 3 



Qo(t) 



f («?Qo(r) - Qo(r)- 1 ) dr 
Jo 



K-2 



-4^(0 I 2)+ Anj (0 \l) + 6iK 2 lK2 Q2 + {2k iK3 + 4)Q 
-2in 2 Qi - k\ (2 I 0) + 2kI (2 I 0) - (l I 0) . (4.5) 
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The functions appearing in the right-hand side of (4.5) are all quasi-periodic, 
since M(Qo) = M(Qi) = M (Q2) = 0. Therefore, we are allowed to take the mean 
value of each individual term above. The result is 



M(J 4 ) = -4^M(0|2) + 4/t?M(0|l)-KfM(2|0) + 2KfM(2|0)-M(l|0). 



By statements (2) and (6) of Proposition 3T and by ( |3.12|) we have 

M(2|0) = -M(0 I 2) = -M(0 |0 |0) = 0, (4.6) 

M(2|0) = M(l |0) = -M(0 |1) . (4.7) 

Hence, 

M{li) = 2k\M (0 1 1) - M (1 I 0) . (4.8) 

We have the following 



Proposition 4.1 Under M{Q ) = M(Qi) = 0, one has M{ 1 1 ) = 2M( 1 1 ). □ 
Proof. By the definition of Qi, 

M(0 |1) = M(0 |0 |0) (4.9) 

and 

M (1 |0) = M (oo(t) (J Qo{r)- l di\ (j Qa{r)- l dn^\ = -M (0 |0 |0) , 

(4.10) 

where, in the last equality, we have used statement (4) of Proposition Now, 
by statement (5) of the same proposition, we can write 

M (0 j j ) + M (0 j I ) + M (0 1 I ) =0 (4.11) 

(recall that M (0 | 0) = M (0 | 0) = M (0 | 0) = 0). Once again, by statement (4) 
of Proposition M (0 | | ) = M (0 | | ) . Thus, flO!] ) reads 2Af (0 | | ) + 
M (0 j I ) = 0. Taking the complex conjugate, yields 



2M (0 I I ) + M (0 I j ) = 0. 

Finally, using and ( p~10l ), we get M (l | 0) = 2M (0 | 1). ■ 

We have just proven that M(J 4 ) = 2 M (Q 1 1 ) — M (0 1 1)J . Now we 

impose M(Ii) — 0. Of course, this will be the case if M (0 | 1) = but, in the 
situation where M (0 1 1) 7^ this can be achieved by fixing ki as (see (4.8)) 

/ m(o|i) \ 1/2 / m(^) V /2 r , 19 . 
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Thus, K\ is a phase: \k\\ = 1. It will be henceforth assumed that M(0 1 1) ^ 0. 
If M (0 1 1) = 0, ki has to be fixed by M(Z 5 ) = 0. We shall not treat this more 
restrictive case here. 

So far, we have verified the absence of secular terms in the series expansion 
( |2.5| ) for g up to order three in e and we have eliminated them from by making 
a especial choice for the value of K\ (given by (4.12)). At this point we would like 
to proceed recursively by imposing M(2 n ) — 0, for all n > 5. This would give the 
correct values for the constants n p , p > 2, and guarantee the absence of secular 
terms in all v n , n > 5. This recursive procedure was used in Q to eliminate 
the secular terms from g in cases (I) and (II). Here, we still have to determine 
the constants n 2 and k 3 explicitly (not recursively) before running the recursive 
procedure. 



4.2 The Absence of Secular Terms in Fixing n 2 



Let us begin by calculating the integrand X5 which appears in V5. Taking n = 5 in 
(2.E) we get Z5 = 2viV4 + 2v 2 v 3 . We first evaluate V1 V4 ex plicitly and then V2V3. 
For viVi, a lengthy computation (see Appendix A.l ) shows that 



VlVi = -Qk\k 2 9* 2 (0 |2)- 



2kik 2 «K 2 (0 1 1) + iKi(2KiK 3 + k\) Q 2 



(4.13) 



where is the quasi-periodic function defined in (A. 7) and depends only on 
the constant K\. Since we are working under the condition M(Qq) = M(Q±) — 
M(Q,2) = 0, we can drop the symbol 9\ 2 above. From this fact and from relation 
we conclude that M (viv 4 ) = 2k 1 k 2 M (0 1 1 ) + imM(Ai). 
Let us now calculate the second term in I5, namely, W2W3. Another lengthy 
computation (see Appendix |A.2[) gives 



v 2 v 3 



-2iKlm 2 (2 \ 2) + 2i K l9\ 2 (2 |1)-2kJk 2 $H 2 (2 | ) + 2i«? JR 2 (1 |2) 

-2 lKl $H 2 (l |1) + 2k 1 k 2 SH 2 (1 |0)-2/s?/sa9la (0 ] 2 ) 

+2kik 2 <n 2 (0 1 1) + i(2Ki^ + Kf/c 3 )Qa - wsQi + k 2 k 3 Q . (4.14) 



Clearly the right hand side of equation (4.14) is quasi-periodic. Since we are work- 
ing under the conditio n M( Qq) = M{Q\) = M(Q 2 ) = 0, we can drop the symbol 
9*2 above and reorder (4.14) in the form 



v 2 v 3 



i(2Kn4 + kIk 3 )Q 2 - in 3 Qx + k 2 k 3 Q - 2in\ (2 | 2) - 2i/ci (1 | 1) 



+2i K \ [(2 |1) + (1 |2)]-2k?k 2 [(2 |0; 
+2k iK2 [(1 |0) + (0 |1)] . 



(0 I 2)] 



(4.15) 



Above we also used 9t 2 (0 |0) = Q 2 . From ( fl,15| ), from the fact that M(Q ) = 
M(Qi) — M(Q 2 ) — and from statements (2) and (3) of Proposition we see 
immediately that M(v 2 v 3 ) = 0. 
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We are now ready to find the value of k 2 by imposing the condition M(I§) = 
0, which guarantees the absence of secular terms in V5. Since M(2s) = 2M (uittj) + 
2M(v 2 v 3 ) = 4kik 2 M(0 1 1) + 2ik 1 M(Ai), we conclude that 

* 2 = - mA :\. (4.16) 

2 2M (0 1 1 ) V ; 



Note that the right-hand side of (4.16) depends on the previously fixed k%. 



4.3 The Absence of Secular Terms in i> 6 . Fixing k 3 

We still have to find k 3 in order to fix recursively all k„'s for n > 4. k 3 will be 
fixed by eliminating the secular terms from Vq, that is, by imposing M(Iq) = 0. 
First of all, we need to write Xq. Using relation ( |2.8| ) for n = 6 we find that 
Iq = 2v 1^ 5 + 2 f2t>4 + f 2 . Let us calculate 2«ii>5. Another lengthy computation (see 
Appendix A. 3 ) gives 



2v lV5 = A 2 -12i4K 3 yi 2 (0\2)+4K 1 K 3 m 2 (p\l] 
+4iki(k 2 k 3 + kik 4 )Q 2 + 2kik 5 Q , 



(4.17) 



where A 2 is the quasi-periodic function defined in ( A.ll| ) and depends only on the 
constants n\ and k 2 - Since 2v\v$ given above is a sum of quasi-periodic functions 
we can take the mean value of each individual term which appears in the right 
hand side of (4.17) and write 

2M(viV 5 ) = M{A 2 ) + 4K 3 mM (0 1 1) , (4.18) 

where, once again, we have used M(Q ) = M(Q 1 ) = M{Q 2 ) — and the identity 
fj)- 

We will now calculate the second term in X5, namely, 2v 2 i>4. We have (see 
Appendix A. 4) 



2v 2 v 4 



A 3 ~^lK 3 m 2 (2 |0)+4 KlK3 ^ 2 (l |0) 

+2i(K 2 K 4 + 2mn 2 K 3 )Q 2 — 2iKiQi + 2k 2 k 4 Q , (4.19) 



where A 3 is defined in ( A.12| ) and depends only on the constants k,\ and n 2 . 

We can now proceed and take the mean value of 2v 2 v 4 from (4.19). Using 
M(Q ) = M(Qi) = M(Ga) = and (O), we get 



2M{v 2 Vi) = M{A 3 ) + A Ki k 3 M (0 1 1] 



(4.20) 



We are almost through with the calculation of Iq. We still need v 3 . Using 
relations ([0])-([L4|), 

v 3 = q- 1 {2im 2 (0\v 1 v 2 ) + K 3 Q } 

= q- 1 {-24^(0 |2) + 2 Kl « 2 (0 \1) + 2i Ki k 2 Q 2 + k 3 Q q } . (4.21) 
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Hence, 



where 



\Qo+4k iK3 [-Klm 2 (0 | 2 ) + JH 2 (0 \1)+ik 2 Q 2 ] -4A 4 . 



At := Qo 1 [«k?5H 2 (0 |2)-iKi«2(0 \l) + k iK2 Q 2 ] 2 . 

Again, A4 is quasi-periodic and depends only on the already fixed constants K\ 
and k 2 . Therefore, using once more M(Q ) = M(Q 2 ) — and (4.6), we get 



M(vl) = 4k 3 K!M (0 1 1) - 4M(_4 4 ) 



(4.22) 



Finally, imposing M(l 6 ) = 0, i.e., 2M(v 1 v 5 ) + 2 M{v 2 v A ) + M(«g) = 0, we 
obtain from our previous calculations (relations ( 4.18| ), ( 4.20 ) and ( 4.22 )), 



«3 



4M{A 4 ) ~ M(A 3 ) - M(A 2 ) 
4kiM(0 1 1) 



(4.23) 



Note that the right-hand side of (4.23) depends on the previously fixed ki and k 2 . 

4.4 The Absence of Secular Terms in v n , n > 7. Fixing n n _ 3 
Recursively 

So far, we have fixed the constants k%, k 2 and k 3 individually. Now we proceed to 
fix recursively all other k„_3 for all n > 7. We have to impose 



M(l n ) 



M ^2 y p v n~p 
Vp=i / 



(4.24) 



for all n > 7. Condition ( 4.24 ) guarantees the absence of secular terms in all v n , 
n>7. 

The idea now is to use ( 1.24 ) to calculate recursively the constants k„_3, for 
all n > 7, that is, K4, K5, .... Of course, we already have «i, k 2 and K3 and, 
hence, we completely know v\, v 2 and V3. We also know that all functions, from 
V\ to vq, are quasi-periodic. From now on we will work inductively. Thus, it will 
be supposed for each n > 7 that we fixed k±, . . . , n n -i by imposing M{X m ) = 
for all m = 2, . . . , n — 1 and that, as a consequence, all functions v\, . . . , w„_i 
are quasi-periodic. Note that it is indeed true for n = 7. 

By our inductive hypothesis, we are allowed to take the summation out of 
the mean value M in (4.24) and write 



2[M(viV n - 1 ) + M(v 2 v n - 2 ) + M(v 3 v n - 3 )] + ^M(v p v n - P ) = 0, (4.25) 

p=4 
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where, by convention, X™=4 M(v p v n - P ) = for n = 7. Let us introduce now the 
following definition: 



l m {t) := q(t){v m (t) - K m q{t)) , 



(4.26) 



for all m > 4. Note that, by relation (2.8), the functions l m 's above can also be 
written as 



J v p( T ) V ™-p( T ) dT 



(4.27) 



For m < n we are allowed to write 

„t m— 1 



* «w 2 / M r ) u "»-p( r ) - M («p 

v° p=l 



m— 1 



m-p l t i 



(4.28) 



p=i 

im- 1 



since we assumed M(Y^Ti v p v m ^ p ) — M(I m ) — for all m < n, by the inductive 
hypothesis. Hence, l m are quasi-periodic for all m < n, by the inductive hypothesis. 

Let us use the definition given in ( 4. 26] ) and evaluate the first three terms 
which appear in (4.25). Beginning with the first one, we have 



M(v X V n - X ) = M(vi(q 1 l n - X + Kn-lS)) = K\M{l n _ x ), 



(4.29) 



where we have used (2.6) and the fact that M(Qo) — 0. Using (4.26), the second 
term of (4.25) can be evaluated as 



M(v 2 v n - 2 ) = M(v 2 {q 1 ln-2 + n n -2q)) = M(q~ 1 v 2 l n -2), 



(4.30) 



where we have used ( |4.2| ) to express qv 2 and t he fact that M(Qo) = M(Qx) = 
M(Q 2 ) = 0. Finally, for the third term of ( |4.25|) , we have 

M(v3V n -3) = M (v3(q~ l l n -3 + K n -3Q)) = M (q'^^Vslns) + K n - 3 M (qv 3 ) . 

(4.31) 



The product qv 3 can be obtained from (4.21), from which we conclude that M(qv 3 ) 
2kxM (0 I 1). Inserting this into j4.3l| ), gives 



M(v 3 v n - 3 ) = M(q- 1 v 3 l n - 3 ) + 2KxK n - a M(0 |1) 



(4.32) 
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The substitution of fll2S| ), ( ^.30Q and ( f4.32| ) into ( |4.25| ), gives us 



2 < Kl A/(Z„_ 1 ) + A/( g - 1 w 2 Z„_ 2 ) + M( (7 - 1 t>3Z, i _3) + 2 K „_3 Kl Af (0 |1) 



(ii) 



(Hi) 



(iv) 



^M(v p « n _ p ) = 0. (4.33) 



p=4 



(v) 

Before we proceed, let us make some comments on our strategy. Equation 



(4.33) is a direct consequence of (4.24) and, thus, is being imposed for each n > 7, 
leading to the values of m, K5 and so on. By our induction hypothesis, we have fixed 
the constants K\, . . . , K n _4 and, hence, we completely know V\, . . . , tVi-4- For this 



reason the terms (hi) and (v) are known by assumption (by (4.28), the evaluation 
of Z„_3 requires the knowledge of V\,... , tVi-4). Our aim is to use (4.33) as a 
condition to fix k„_3 and we, therefore, have to isolate the dependence of (4.33) 
on K n -3- The function l n -2, in ter m (ii), depends implicitly on i>„-3 and, hence, 
on n n -3 (see, again, relation ( 4.28 )). The term (i), however, depends implicitly on 
k„_2 and K n -3- This dependence on k„_2 could be a problem, since we are still 
working to fix k„_3. Nevertheless, as will be shown, the conditions M(Qi) = 
and M(Q,2) — fortunately eliminate k„_2 from the hnal expression, and we will 
be led to a condition expressing k„_3 in terms of known quantities. 

Let us now compute terms (i) and (ii). After a long computation, found in 
Appendix A. 5 , we get 

KxM (l n -l) = 2KiK n . 3 M (0 I 1 ) + TlW , 



(4.34) 



where IZn^ i s defi ned in ( A. 18 ) and depends on constants K\, ... , n n -i only. For 
term (ii) of (4.33) we get, after another long computation presented in Appendix 



A.6 



M{q- 1 v 2 l n _ 2 ) 



-2 Kl K n -3M{0 \ i)+n 



(2) 



(4.35) 



where TZ^ is defined in (|A.19|) . We again stress that Teh depends on K\ , 
only. 

We ar c now read y to g ive the precise value of k„ _3 in order to satisfy ( 4.24 ) 



,(2) 



Collecting ( 4.34 ) and ( 4.35 ) and inserting them into ( 4.33 ), we obtain 
-1 



4 Kl M(0 1 1) 



M(v p v r , 



2fcW + 27^ 2) + 2M{q- 1 v 3 l n -s) 



for all n > 7. Note that izffl and can be completely computed in any order 



of recursion (see equations ( A. 18 ) and ( A. 19 )). The f uncti on l n -3 can also be 
computed in any order of recursion if one uses relation (4.28). 
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Summarising our conclusions, under conditions (IIIo) and M (0 | 1) ^ 0, i.e., 
under condition (III) of Theorem |2.2| , and with the constants K n recursively chosen 
as 

1/2 

(4.36) 



Ki = 



K2 



M (0 | 1 ) 



\ M(0 1 1) / 

2M(0 1 1) ' 
-1 

4kiM(0 1 1) 



«3 



4Af (AQ - M(^ 3 ) - M(^ 2 ) 
4kiM(0 I 1) 



(4.37) 



' n— 4 



^ M (u p v n _ p ) + 27^ 1} + + 2Af ((TVn-a) \ , 

.p=4 J 

(4.38) 



for n > 7, all secular terms are eliminated from the formal solution (2.5) of the 
generalised Riccati equation ( |2.1[ ) . Notice that hypothesis M (0 1 1 ) 7^ is the 
only additional restriction needed to (4.36)-(4.38). The proof of Theorem is 
thus complete. ■ 



5 The Case of the Monochromatic Field 

We illustrate our method and our results considering the simplest case of monochro- 
matic interactions (ac-dc field) 

f(t) = F + <pcos(ut), (5.1) 

which are of particular importance in physical applications. We want to show 
that with conditions (I)- (III) we obtain with our method convergent perturbative 
solutions of this problems for all parameters Fq and <p, except perhaps for some 
spurious situations. For ( ]5.1| ) one has 

Go(t) - E J »to) el(n+X2M ' 

where J n is the Bessel function of first kind and order n and where we defined 
Xi ■= 2ip/cj, X2 '■= 2F /cj. Hence, condition (I) (treated in detail in ||) holds 
provided X2 = —m, with m int eger, and provided xi is n °t a zer0 of the Bessel 
function J m . Here, we have by ( 2.14| ), 



= -— + eJ m ( X i) + 0(e 2 ) . 

See also the discussion in [Q. Let us consider the complementary situations. 

i. Consider the case where X2 = —m, a non-zero integer, and xi is a zero of 
J m . One has M(Qq) — and we have to look first at M(Qi). We get, 

Af(Gi) = X Jk+m ^ 2 



~i iku> 
feez* 
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For integer m one has 



E 



d 

-2—J„{x) 
ov 



+ irY m (x) 



(5.2) 



where ^Y m {x) = j^r,J v {x)- (-l) m g^ J- V (x)) \ v=m , Y m being Bessel functions 
of second kind. Identity fl5.2|) i s apparently unmentioned in the literature, and we 
sketch its proof in Appendix |D|. Since xi is bound to be a zero of J m , one concludes 
from this identity that M(Qi) = in this case. A direct computation shows that 



M (0 1 1) = -^T m ( X i), where T m (x) := - ^ 



k,p^0 



Jra+p{x)Jrn+p—k {x)Jra+k{x) 

kp 



Numerical calculations indicate that T m {x) does not vanish at the zeros of J m . We 
refrain from giving an analytic proof of this fact. We conclude that condition (III) 
holds in case i, except, perhaps, for spurious zeros of J m for which % n even tually 
vanishes, and whose existence could not be ruled out numerically. By (2.15), 



raw 



2c 



T m ( Xl )+0(e 4 



ii. Consider the case where \2 is non-integer (see also the discussion in 
One has M(Qo) = and we have to look at M(Qi). We get, 



M(Qi 



i 

UJ\2 



^o(xi) 2 + 2^E v 2 

fc=l X2 



■hixif 



k 2 



Generally the r.h.s is non-zero and we have condition (II). Hence, by ( p. 15 ) 



n = F 



^X2 



Mxif + txlY, 



■Mxif 
^ xi-k 2 

fe=l A2 



+ 0(e 3 ). 



(5.3) 



Jfc(Xi) 

U2~ T 



Note, however, that on each interval X2 € (fc, fc + 1), fc = 1, 2, . . . , the terms 
— r ^^(^Y^ vary continuously from +oo to — oo. Hence, there is on each 
interval (fc, k + 1), k = 0, 1, 2, . . . , a special value X2 °f X2 (depending on xi) 
for which M(Qi) = 0, and we would be out of case (II). But when 2\2 is non- 
integer, one has M(Q^) = 0. Hence, except for the very unlikely case where 
is an integer, we would be out of condition (III) as well, and Q = F + 0(e 4 ). 
This special case may not be very interesting, because X2 has to be chosen with 
precision. That could be difficult to fix it in some experimental setting. 

Another special situation would occur when \2 is chosen to satisfy Fq — 
ie 2 M(Qi) — 0. By the argument above, this is possible, but X2 will depend on 
e. It is therefore unclear if £7 will be just 0(e 4 ) or "small" (eventually leading to 
an even stronger dynamical localisation than we have in case (III)). It is not even 



2G 



clear that we will be in a situation where our series converge, and we left this other 
special situation without more comments. 

It is interesting to compare the expressions for the secular frequency Q in 
the three situations above (for Fq ^ 0) with the situation where ip = 0, where the 

secular frequency f^o is Oq := Fq./i + (J-) 2 = F + £r + 0(e 3 ). This reveals the 



effect of the ac-field ipcos(ojt) on the secular frequency. Taking xi in (5.3) we 
recover Qq. 

iii. Consider the case where %2 = 0, i.e., Fq = 0, and xi = x a, the a-th zero 
zero of the Bessel function Jo on the positive real axis. This case is interesting in 
connection with the issue of dynamical localisation, as discussed by many authors 
(p^j. For more references, see §, |). Here M(Q ) = and, 

M(Ql) = jy(Jm(Xl) 2 -J-m(Xl) 2 \ = „ 
m=l v ' 

since Jk(x) = (— l) k J-k(x). Thus, condition (II) does not apply and we have to 
look at M (0 1 1). We obtain M (0 1 1) = -4-T (x a ), with 



Jn (-^a) J n—m (^a) Jm (*^a) 

nm 



We conclude that condition (III) will be valid, except perhaps for spurious zeros of 
Jo for which T (x a ) eventually vanishes. Numerical computations, though, indicate 
that such zeros may not exist. We abstain from giving an analytical proof that 
T (x a ) 7^ for all zeros of the Jo function. Figure |l] shows the values of |T(ai a )| 
calculated numerically on the first fifteen zeros x a of Jo. We note an exponential 
decay of |T(x a )| as x a increases, suggesting that \T (x a ) \ — * only for a — > oo. 

We conclude that condition (III) is suitable for studying the monochromatic 
field when xi li es over the "resonant" points x a , leading, together with condi- 



tion (I), to a complete solution of (5.1) except, perhaps, for some rather spurious 
situations. 



Note finally that in this case iii we have fl = 0(e 3 ) (see ( p. 16 )). In fact 



the first contribution to 17 is (x a ). This weak dependence on e implies long 

transition times for certain probability amplitudes (see below), an issue known as 
dynamical localisation (see 0). 

In order to test our algorithm and to extract more information from our 
solutions, we wrote a computer program to compute the matrix elements of the 
propagator U(t) given in (^) for the case iii described above, where one has ap- 
proximate dynamical localisation. The results are excellent and are briefly reported 
in the next section. 

Remark 5.1 We obtained for case iii, f2 = (x a ) + 0(e 4 ). Following another 
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0.0 10.0 20.0 30.0 40.0 50.0 

Figure 1: The filled circles represent the values of |T (x a )\ as a function of the a-th 
zero x a of Jq in R + . We infer an exponential decay of \M (0 1 1 ) | as a increases. 



path, we obtained in 



t J2n 1 + l(Xa)J2n 2 + l(Xa)J-2(n 1 +n 2 + l)(Xa) , n( 4^ 

^ (2n 1 + l)(2n 2 + l)] + U *>- 

Til, n 2 GZ 

Thus, we must have 

Jn (*£a) Jn—m ip^a) Jm (*^a) 



E 



J2ni+l{ x a) J2n 2 + l{ x a) J -2{ni-\-n 2 -\-l) Oa) 



E 



(2711 + 1) (2712 + 1) 

Curiously, such an identity is not easy to prove by direct means. It can be estab- 
lished, however, with the use of statement (6) of Proposition 3.1, namely, from 
M (0 | | ) = 0. Actually, such an identity is not specific for Bessel functions, 
but is valid for sequences satisfying certain properties. 



6 Numerical Results 



In this final section, we will briefly show some numerical results that were obtained 
using condition (III) to study the monochromatic interaction over the resonant 
points (see the discussion in the last section). O ur g oal here is to show how one 



can numerically use the results of Theorems 2.1, 2.2 and 2.2 to produce results of 
physical interest. 

The com putation of the wave function <&(<) associated with the Schr odin ger 
equation (L3) can be performed by means of the propagator U(t), given in (2.2), in 
terms of the solution of the generalised Riccati equation (2.1) given in Theorem 2.2. 
Due to the uniform convergence of g(t) for all t e M (see Theorem UM, one can use 
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the proper expressions given in to compute the elements of U(t) via uniformly 
convergent power series in e. Although this task may not be so trivial, it can be 
implemented numerically with great success, as we are about to show now. 

The first step toward the computation of U(t) is to evaluate the Fourier 
coefficients of g{t) using the recursive relations given in (B.4)-( |B~2] ). Of course, we 
need first to compute the constants n n . To this end, we have to use expressions 
(4.36) and ( 4.38| ). Those equations involves the mean value of some previously 
defined functions, which can be obtained trough the zero order terms of their 
Fourier decomposition. Equation (3.4) of Proposition 3.4 can be used to this end. 
Once Vm^ 1 are all computed, the m-th Fourier coefficient of order n of g(t) can 
be obtained via the convolution Gm' = EpexQra-p^p"' ( see equation (|2.5|)). 
As a second and final step towards the computation of U(t), we have to use the 
expressions given in Basically, these expressions convert the integral form of 
the elements of the propagator U(t) (see Theorem |2.4[ ) into convergent power series 
in e using the Fourier coefficients Gm' of g{t). They also lead to the expression of 
U(t) in terms of its Floquet form (see Theorem [Lf]). A more detailed description 
of our numerical study will be postponed for a future publication ]13| ], For now, 
let us briefly show some of our results. 

We have followed steps one and two above to numerically compute the prop- 
agator U(t) associated with (4.3) when f{t) = (pcos(uit) in the situation where 
M (Qo) = 0, that is, over the "resonant" points defined by the condition (p = 
ux a /2, x a being the a-th positive zero of Jo- 

Let = (o) an( l ^- = (l) ^ e t w0 orthogonal states of a system described 
by (4.3). The probability for the transition from the initial state $+ to the final 
state at time * is P(t) := | ($+, U(t)$-) | 2 = \U 12 {t)\ 2 . We computed P(t) 
numerically using our expansions. To estimate the accuracy of our calculations, 
we tested the unitarity of the time evolution and considered the quantity N(t) :— 
|^n(i)| 2 + |C^i2(i)| 2 , which should be identically equal to 4 for unitary U(t). 

In Figure 0(a), we show P(t) for ip = uxi/2, x\ being the first positive zero 
of J . We took uj = 40.0, e = 0.4 and worked with a sixth order expansion in e. 
The time interval considered corresponds to 4.6 x 40 6 times the basic cycle 2tt/u> 
of /. We recall that in this case we have condition (III) and the secular frequency, 
which dominates the quasi-pe riod ic evolution of the system, is of order e 3 /w 2 (see 
H or the discussion in Section 2.1). This explains the long time (~ 40 6 of the basic 
cycle of /) needed for the system to transit from $ + to «3>_. 

In Figure ||(b), we show the quantity N(t) — 4 as a function of time obtained 
from the same calculations leading to Figure |](a). By looking at the deviations 
of N{t) from 1, we can infer that our perturbative solution produces errors of 
the order of only 0.3%. This excellent numerical precision for long times is a 
consequence of the uniform convergence in time of our expansions, i.e., of the 
elimination of the secular terms. 

In Figure H we took lo = 1.0, e = 0.04 and (pi = uxi/2. From Figure |3](b) we 
infer errors of the order of only 0.05%. The time interval considered corresponds to 
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Figure 2: (a) The transition probability P(t) as function of time. Here e = 0.1, 
lo = 10 and 2(pi/u> — x\, the first zero of Jo. (b) The quantity N(t) — 1 that 
measures deviation from unitarity. 

1.6 x 10 6 times the basic cycle 2it/uj of /. In Figure || we took ui — 10.0, e = 0.01 
and again ipi = uix\/2. From Figure ||(b) we infer errors of the order of only 
0.003%. The time interval considered corresponds to 1.6 x 10 9 times the basic 
cycle 2it/uj of /. 

The numerical computations above involved the evaluation of the Fourier 
expansions of several functions, like q, i> n 's etc. We typically computed about 20 
terms of that expansions, but since the coefficients decay very fast, even less terms 
are needed. More details about the numerical computations will be presented in a 
future publication 




0.0 2.5 5.0 7.5 10.0 0.00 2.5 5.0 7.5 10.0 

Time [ x 10 6 ] Time [ x 10 6 ] 



Figure 3: (a) The transition probability P(t) as function of time. Here e = 0.01, 
lu = 1.0 and 2ip\jio = x\, the first zero of Jo. (b) The quantity N(t) — 1 that 
measures deviation from unitarity. 
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0.0 2.5 5.0 7.5 

Time [ X 10 8 ] 



10.0 



0.00 



2.5 5.0 7.5 10.0 

Time [ x 10 8 ] 



Figure 4: (a) The transition probability P(t) as function of time. Here e = 0.01, 
lo = 10.0 and 2<p\/u) — X\, the first zero of Jo. (b) The quantity N{t) — 1 that 
measures deviation from unitarity. 



Appendices 



A Some Special Relations 

This appendix presents the proofs of some relations used in Section |]. Since they 
involve a somewhat large amount of algebraic manipulations we prefer to separate 
them from the main text. In the following we will (often without explicit men 
tion) make repeated use of the propositions and corollaries proven in Sections 3.1 
and EDI. 



A.l Obtaining Relation ( |4A3| ) 

Since we have Mil?,) = and since with the choi c e of n\ in ( 4.12| ) we imposed 
M(Ii) — 0, we have by the recursive relations (4.1)-(4.3) 



V!V 3 = 2ini 9\ 2 (0 | v x v 2 ) + kik 3 Q , (A.l) 
V1V4 — 2in\ EH2 (0 I ) + iki <K 2 (0 \v 2 2 ) +K lKi Qo. (A.2) 



Inserting (4.4) into the r.h.s. of (A.l) we get 
Vl v 3 = -2 K j d\ 2 {0 \2) + 2 K f m 2 {0 \l) + 2i K 2 lK2 m 2 (0 \0) + kik 3 Q . (A.3) 



Inserting this into the r.h.s. of (A.2) gives 

ViV 4 = 



-4i/sf (0 I I 2 ) + Ainl m 3 (0 1 1 1 ) - A k \k 2 «R 3 (0 | | ) 
+2iK\K 3 y\ 2 (0 |o) + i Kl m 2 (0 I v\ ) + K1K4Q0 . (A.4) 
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Let us now compute v\. Since v 2 = k 2 9 + icT 1 {k\Q,2 — Qi), we have 

V 2 = h\Qq + 2«KjK 2 Q 2 — 2ik 2 Qi + s , 

where s := — Qq 1 (ki<2 2 — Si) 2 - N ote that s contains k,\ alone. 
Now, we insert this into (A.4) and get 



(A.5) 



Vl v 4 = -Un\ ^3(° 1° |2)+4i K J 5H 3 ( |o 1 1 ) - 6k?k 2 ^2 (0 |2) 
+2i«?/s 3 ^2 (0 |0) + i«i«| 3t 2 (0 |0) + 2 Ki k 2 9t 2 (0 |1) 
+iKi9l 2 (0 |s) + kik 4 Qo, (A.6) 



where we also used the fact that 9t 3 (0 | | ) = 9t 2 (0 | 2). 

Collecting in (|A.6|) the terms depending only on Ki, we define 



(A.7) 



A x := -A4 9t 3 (0 I I 2 ) + 4k 2 m 3 (0 | 1 1 ) + <K 2 (0 | s) 
and rewrite ( |A.6| ) as 

v lV4 , = -6k?k 2 912 (0 \2)+iK 1 (2KiK3+i4) Q 2 + 2k 1 k2 m 2 (0 | l) + K 1 K 4i Qo+iKlAx, 



which is expression ( 4.13 ), as desired. From definition ( A.7 ), it is evident that Ai 
is quasi-periodic, since it is a sum of quasi-periodic functions. 

A. 2 Obtaining Relation (|4.14j) 

From relations 



3)-dyD 

V2V3 = Qo 1 (iKlQ 2 - iQi + k 2 Qo) (2im 2 (0\v 1 v2) + n 3 Q ) . (A.8) 
Now, by (Q, we have Qq x Qj 5R 2 (0 | v lV2 ) = 9t 2 (i | fi« 2 ). Hence, ([Aj) becomes 
t> 2 ?; 3 = -2k?$K 2 (2 |w 1 w 2 ) + 2fR 2 (l |u lU2 ) + 2iK 2 fR 2 (0 |viw 2 ) 



-«3 (iKiQ.2 - iQi + k 2 Qo) 



(A.9) 



Inserting (4.4) into the r.h.s. of (|A.9|), it becomes 



v 2 v 3 = -2i«$SRa(2 |2) + 2i«?JRa(2 1 1) - 2k?k 2 JH 2 (2 1 0) + 2i«? 9ta (1 |2) 
-2iKi«Ra(l |1) + 2k 1 K2«K 2 (1 | 0) - 2k?k; 2 9t 2 (0 |2) 
+2«ik 2 JH 2 (0 I 1 ) + i(2«iK 2 + k 2 k 3 )Q 2 - in 3 Qi + k 2 k 3 Qo . 



This is (4.14), as desired. 
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A. 3 Obtaining Relation (|4.17|) 

According to equations (4.1 )— (|4.3D , we have 

2vxv 5 = Hki^ 2 (0 \vxVi + v 2 v 3 ) +2/tiK 5 Q - (A. 10) 

By ( §3f ) and 

+ = -2mf£H 2 (2 |2) + 2i/t?£H 2 (2 1 1 ) - 2k?k 2 % (2 |0) 

+2ire?9t 2 (l |2)-2iK 1 $H 2 (l 1 1 ) + 2 Ki k 2 1R 2 (1 |0) 
-8k?K2^2 (0 |2)+4k iK2 « 2 (0 \l)+3i(K 2 1 K 3 + K 1 K 2 2 )Q 2 
— «k 3 Qi + (k 2 k 3 + KiK 4 )<2 + i/«i-4i • 
Inserting this into ( A. 10] ), gives 



2v lVa = 8re?<H 3 (0 | 2 | 2 ) — 8i/cf 9t 3 (0 |2 | 1 ) - 8mf K 2 «3 (0 |2 |0) 
-8 K f9t 3 (0 1 1 | 2 ) + 8«? SR3 (0 1 1 1 1 ) + 8ik\k 2 m 3 (0 1 1 | ) 
-32mf«2 % (0 | 1 2 ) + 16ik?k 2 fR 3 (0 | 1 1 ) 

-12(^ K3 + K 2 K 2^ 2 (q | 2 ) + 4K!K 3 0^2 (0 1 1 ) 

+4iK 1 (K 2 K 3 + KiK 4 )fR 2 (0 |0)-4K?9t 2 (0 | ^i) + 2k 1 k 5 Q ■ 

Let us now isolate all terms depending only on k± and/or n 2 in one single term, 
which we call A 2 : 



A 2 



8/c?m 3 (0 | 2 | 2 ) - 8i K fd\ 3 (0 | 2 1 1 ) - 8i^K 2 5K 3 (0 | 2 | ) 
-8/4 JR 3 (0 1 1 | 2 ) + 8/c? SR 3 (0 1 1 1 1 ) + 8i/e?K 2 <K 3 (0 1 1 | ) 
+16iK 2 lK2 £Ra (0 | 1 1 ) — 32iKjn 2 W 3 (0 | 1 2 ) 
-12k?k3 £H 2 (o 1 2 ) - «n 2 (o | Ai ) 



(A.ll) 



and write 



2vxv 5 



A 2 - 12k?k 3 Ui 2 (0 I 2 ) + 4 Ki k 3 JR 2 (0 1 1 ) 
+ 4jki(k 2 k 3 + kik 4 )1H 2 (0 | 0) + 2kik 5 Q . 



This is ( 4.17 ), as desired. Note that A 2 is clearly quasi-periodic, being a sum of 
quasi-periodic functions. 



A. 4 Obtaining Relation (|4.19|) 

Using relations @)-(|J), and fT|) (which implies Q a 1 Q l $K 2 (0 | a) = Dl 2 (i \ a)) 
we have 

2v 2 v± = 2Q-Q 1 {ik\Q 2 - iQ l + n 2 Q ) [2i9t 2 (0 \ Vl v 3 )+im 2 (0 |v 2 )+k 4 Q ] 
= -4k?«R 2 (2 |wiw 3 )+4$K 2 (1 | i>iv 3 ) + 4ik 2 d\ 2 (0 |«i« 3 ) 



2k? % (2 



-29t 2 (l 



2iK 2 K 2 (0 



-2k 4 {in\Q 2 - iQi + k 2 Q ) 
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Inserting above the expression for V1V3 given in ( |A.3| ), we get 

2v 2 V4 = 8ref«R 3 (2 |0 |2)-8Kf«n 3 (2 |0 |1)- 8^2^3(2 |0 |0) 
-4k?k 3 5H 2 (2 |0)-8/c?9t 3 (l |0 |2)+8k?$K 3 (1 |0 |1) 
+8m?K 2 5H 3 (l |0 |0)+4 Ki k 3 <K 2 (1 1 ) - 8*«f k 2 % (0 |0 |2) 
+8ire?/c 2 9*9 (0 1 1 1 ) - 8k i k 2 ^3 (0 I I ) + Ai Kl K 2 n 3 SH 2 (0 1 0) 
-2k? % (2 1 ) + 2 % (1 I «|) + 2iK 2 9* 2 (0 | «|) 
+2k 4 (ik?Q 2 - iQi + k 2 Qo) • 

We now isolate all terms depending only on n\ and/or k 2 in one single term, 
which we call A3: 

A 3 := 8«f9l 3 (2 |0 |2)-8«f9l 3 (2 |0 |l)-8m^ 2 %(2 |0 |0) 

-8k? 9*3 (1 1 I 2 ) + 8k? 9*3 (1 1 1 1 ) + 8tK i K ^ 9t 3 (1 | | ) 
-8^k 2 9t 3 (0 I I 2 ) + %iK 2 x K 2 9t 3 (0 I 1 1 ) - %k\k\ fR 3 (0 | | ) 
~2 K lm 2 (2 \ v%) +2SK 2 (1 I v\ ) + 2iK 2 m 2 (0 \ vf) (A.12) 

(note here that, by ( |A.5[ ), v\ depends only on K\ and k 2 ). We can now write 

2v 2 v 4 = ^3-4k?k 3 JR2(2 |0)+4 Ki k 3 ^ 2 (1 |0) 

+2i(nln4 + 2kik 2 k 3 )Q 2 — 2zk 4 Qi + 2k 2 k 4 Q • 



This is the desired relation (4.19). Again, note that A3 is quasi-periodic, for it is 



a sum of quasi-periodic functions. 

A. 5 Obtaining Relation ( ggg ) 

Let us explicitly write fti£ n _i in term (i) of ( [4.33 ) using equation (4.27). We have, 



n-2 

(t) - t/si ^ 9l 2 (0 I ) =: 

p=i 



where 



A(t) := 2i «i 9l 2 (0 | viv n - 2 ) , B(t) := zki 9i 2 (0 | ). (A.13) 

p=2 

The above expressions for A(t) and B(t) will now be worked indiv idu ally. L et u s 



start with A(t). By the inductive hypothesis, we are allowed to use (LI) and (4.3) 
We write 

n-3 \ 
VpVn—2—p ) T K n _ 2 Qo 

p=i / 



A(t) = 2 i k{ m 2 ^0 

n— 3 

= -2 K ?^ 9^3(0 |0 |vp«„- 2 -p) + 2i«2 Kn _ 2 fK 2 ( |o) . 



0=1 
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Note that A{t) is implicitly dependent on k„_ 3 , namely through v n _ 3 . To 
make this dependence explicit we have to split the sum containing u n — 3 and write 
u„_3 with the use of (4.3): 

ri— 4 

A(t) = -2/^5^ 913(0 |0 |w p « n _2- P )-4/s?9la(0|0|t;i» n _3) 

p=2 

+2iKlK n - 2 m. 2 (o 1 0) 

n— 4 n— 4 

= -2 «3 (0 I I v p v n - 2 ^ p ) - 4t 4 J2 ^ 4 (° I I I V«-3- P ) 

p=2 p=l 



-4K? Kn _3fH 3 (0 |0 \ 0) + 2iK 2 1 K„- 2 ^2 (0 |0) 



(A.14) 



We will now work on B(t), equation (A. 13). Using (4.1) and (4.3), we have 

n— 4 

B{t) = 2i Kl 9* 2 (0 I v 2 v n „ 3 )+i Kl J2 ^2 (0 I VpVn-t-p ) . (A.15) 

p=3 



Next, we have to compute separately v 2 v n _ 3 . Using once more (4.2), (4.3) and 
(O) (which implies 2(5" 1 (0 | a) = W. 2 (i \ a)), we get 



V 2 V n -3 



Qo 1 (iK(Q 2 - iQl + K 2 Qq) ^2 (0 |»p«n-3-p) t + K„_3Go(t)| 

n— 4 n— 4 

-K\ ^ ^2 (2 I Up w «-3-p ) + ^ 9t2 (1 I 



p=i 

n-4 



p=l 



+««2 ^ 9^2 (0 I UpW n _ 3 -p) + K n -3 {i^iQl - iQl + K 2 Qq) 



0=1 



This expression for v 2 v n — 3 has to be introduced into the first term of ( A.15| ). The 
result is 



n— 4 



fl(t) - 2i/«i ^ < -K? SH 3 (0 I 2 I u p «„_3_p ) + 9t 3 (0 1 1 | 



p=i 



n— 4 



+*«a £H 3 (0 j I ) [ + i«! J] 9t 2 (0 I 

J P=3 

+2m 1 / t „_ 3 [*«i SRa (0 I 2 ) - * 9t 2 (0 1 1 ) + k 2 Ui 2 (0 I )] . (A. 16) 

Since both A and B are quasi-periodic, we can now compute K\M (l n —i) = 
M(A) + M(B). Using JO^) , ( |A~16| ), ([□}) and the already proven fact that 
M (0 |0) = 0, we get 



KiM(l n —i) = 2/si/s n _ 8 M(0|l) + 724 1) , 



(A.17) 
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where 



n— 4 i 



:= 2i Kl J2\ -2k?M(K 4 (0 I I I v p v n . 3 . p )) + M(9* 3 (0 1 1 | u P u„-3- P )) 
P =i { 

— k 1 M( < R 3 (0 | 2 | )) + i/t 2 M( 9t 3 (0 | | w p w n _3_p )) > 



-2k?^M(9* 3 (0 |0 |w P v„-2-p))+iKi 5Z M (^2 (0 |WpW„_l_p)). 

p=2 p=3 

(A.18) 



This is the desired relation (4.34). By inspection, one verifies that IZn^ depends 
on the constants «i, ... , k„_4 only. Notice that the constant K n _2 disappeared 
completely when we took the mean value of A(t) + B(t), due to the crucial fact 
that M(Q 2 ) — 0. This is very important, otherwise we would have in (A. 17) an 
equation for two unknowns k„_ 3 and K ra _2- 

A.6 Obtaining Relation (gj5g) 



The mai n point is to make the K n _3 dependence of £ n _ 2 (i) explicit. Using ( 4.26 
and (4.3) for v n s(t), we can write 



n— 3 



n— 3 



q 1 v 2 l r , 



-2 = n 



~ 1 V 2 ^ 9* 2 (0 | V p V n ~2~p ) = VpV n —2—p ) 



0=1 



p=l 



n— 4 



2? 9* 2 (9^2 I Wl«n-3 ) + i ^ ^ 2 ( qV2 I u P w n-2-p ) 
p=2 

n— 4 



2lKi 9*2 <7«2 



9*2 (0 I WpU„-3-p ) + i^n-iQa 



p=l 



+i 22 9*2 (qV2 I V p V n ^2-p ) 
p=2 

n— 4 n— 4 

= -2ki ^3 O2 I I ) + j^9* 2 {qv 2 \v p V n -2-p) 
p=l p=2 

+2i/tiK„_ 3 9*2 (<?«2 |0) . 



According to (f4.2|), 

9*2 (qw 2 |0) = i«?9t 2 (2 |0)-i9t 2 (l |0) + k 2 9*2(0 |0) 
and, hence, M(9* 2 (qv 2 | 0)) = -iM (1 | 0) = iM (0 j 1). Therefore, 
Miq-Wn^) = -2kik„_ 3 M(0 |l) + ^ 2) , 
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where 



:= -2m ^ M(m 3 (qv 2 | | v p v n - 3 -p )) + M(m 2 (qv 2 \ v p V n - 2 - P )) . 

p=l p=2 

(A.19) 

I k (2) 

This is the desired equation ( [4.35[ ). By inspection, one verifies that 1Z n depends 
on the constants «i, ... , K n _4 only. 



B Proof of Convergence of the e Expansion for 
Periodic / 



Here we will present the proof of Theorem 2.3, i.e., the proof of convergence of 
the e expansion of (2.5) for periodic /. It follows the ideas of but technical 
adaptations are necessary. For the sake of simplification we shall consider here 
only the case where Fo — M{f) = 0. The general case Fq ^ can be treated 
following the lines described in detail in . 

(21 r- 1 

In terms of the Fourier coefficients Q m and Q m , appearing in fl2.11|) , of the 
Fourier coefficients Vm of ( 2.12 ) and of the constants k„, relations (4.1)-(4.3), 
become 



* m 

y(n) 



riiuj 



E 



Qm— (711+712) 

(ni + n 2 )u 



y(p)y( n -p) 

/ j "l «2 
\p=l / 



+ K 2 Q m , (B.l) 

, for n > 3. (B.2) 



Of course, due to the choices of the constants K n described before, no secular terms 
appear. 

By ( 2.13| ) and by an inductive argument, we will prove the following state- 
ment: for all p € N and all m £ Z there are constants K p > such that 



y(p) 



-X\m\ 



To show this, let us first recall the following result, proven in 



(B.3) 



Lemma B.l For x > and m 6 Z define 

B(m) = B(m, x) '•= ^ 



-X(\m-n\ + \n\) 
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Then 



has 



B(m) < B - 



-X\m\ 



<m> 2 ' 

for some constant Bq = -Bq(x) > an d f or a ^ m€Z. 



□ 



From (B.l) and ( 2.13 ), we have 

e -x\m\ 



* m 
V (2) 



< Q- 



< 



2Q 2 



E 



s-X(l"»-"l| + l"l|) 



w i «m — ni» 2 <§;«i3> 2 |m| 



l«2ie 



e -xM 
<m> 2 



where we used = 1. Now 

e -x(l"i-™il+l«il) i 



E 



<m - ni> 2 <?ii> 2 m 



< 



E 



e -x(l™-«i|+|ni|) 
<m - ni> 2 <?ii> 5 



< B - 



-X\m\ 



where the last inequality comes from Lemma [B.l|. Hence, we can write 



V (i) 



<m> 2 
for all m € Z, by choosing 



and 



v^( 2 ) 



-Xl"»| 



2 ' 



and 



X 2 s= ^ + | K2 |Q. 



To proceed, let us assume (B.3) for all p = 1, . . . , n— 1. By (B.2) and (2.13), 
we have 



< 



Q 



( 



E 



-x(l*n— ni— nal+lml+lnal) 



1 



■ » ltn2 ez 



<m — Tii — n 2 > 2 <C?ii> 2 <Cn 2 > 2 |th + n 2 | 
x e -xM 



\p=l 



for n > 3. 



Again, applying twice Lemma B.l, 

p -x(|m-ni-n2|+|ni|+|«2|) 

E 

»1.»2« 



<Cm — ni — n2> 2 <Jii> 2 ^n2> 2 \n\ + n 2 | 



< 



E 



ni, niSL'L 



e — x(|m-m— nal + lml + lnal) 
<m - m - n 2 > 2 <ni> 2 <n 2 > 2 



< (S 



<m> 2 
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Therefore, 



for all m E Z, by choosing 



-X\rn\ 
<?71> 2 



Kr, 



(B.4) 



where k° is some suitably chosen upper bound for \n n \. We now turn our attention 
t o \n n | , for which we have to find estimates using again the inductive hypothesis 
(B.3) for all p = 1, . . . , n — 1. Th e constants Ki, K2 and K3 are are fixed by ( 4.3C| ) 
and K n , n > 4, are given by (4.38), from which we get 



k < 



l 



4|M(0 1 1)| 



n-l 

J2\M{v p v n+3 „ p )\+2 

p— 4 



ft 



(i) 

n+3 
T 2 



+2 



ft 



(2) 
n+3 

T 3 



+2|M(<rVn)| 



(B.5) 



for n > 4. We have to bound each of the terms Tj indicated above. Let us start 
with T\. 



Bound for 7\. By @, M (w p w„+3-p) = £ roeZ V#Vi^ +3 - p) . Hence, by 



the inductive hypothesis (B.3), assumed for p = 1, ... , n — 1, we have 



where 771 := 



n— 1 n— 1 

n+3— p)\ 5~ Tjl ^ ] KpK n +3—p, 

,-2x|m| 



(B.6) 



m£2 



Bound for T2. Expression ( A.18| ) involves sums of the mean value of func- 
tions like 9\k (ai | • • • | Uk )■ Let us prove a general statement about such functions. 



Proposition B.2 For k > 2, Zef ai, . . . 

frequency u>, and such that ai(t) = J2 m ei 
dents (m) satisfy 



A {l \m) 



Ofc be periodic functions with the same 



A^(m)t 



< OLl- 



-X\™.\ 

«m> 2 



where the Fourier coef 



for all m E Z and all I = 1, . . . , k, where a; > and x > 0. 
is a positive constant flk such that the Fourier coefficients D\k {a>i \ 



(B.7) 

Then, there 
■■ \ ak) (m), 
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m € Z, 0/ SHfc (ai 



I dk ) t are bounded by 



1 9tfc (ai I • • • I afc ) (m)| < /3fc «i • • • a k 



e -x|m| 
<m> 2 



(B. 



□ 



Proof. Let us first consider the case k = 2. The Fourier coefficients of 9^2 («i | o, 2 ) t 
are given by 



where 



K 2 |a 2 )(m) = ]TA (1) (r 

nGZ 

f A( 2 )(m) 



I( 2 )(m) := < 



imuj 



A( 2 )(fc) 



— ' fc 

fcGZ, 



for m^O, 



for m = 



(B.9) 



From (B.7), it follows that 



A? 5 (m) 



< a 2 



D e-xW 
uj <m> 2 



where D := £ meZ %S3F 
one has 



Therefore, from (BJ3), by (B7) and by Lemma 
B D e~x\ m \ 



B.l 



1 9^2 (ai \a 2 )(m)\ < a x a 2 - 



This proves the statement for k — 2. The general case follows from (3.9), by 
induction. ■ 
As a corollary, one sees from ( 2.13 ) that the functions Qq, Qi and Q 2 have 
Fourier coefficients bounded as |Qi(m)| < 7i ^m>2 for some positive 7,. Moreover, 
by the inductive hypothesis (B.3) and by Lemma B.l. the Fourier coefficients 
(v p v q )(m) of product functions like v p (t)v q (t), with p, q = 1, . . . , n — 1, are also 
bounded as 



\(v p v q )(m)\ < B K p K q 



e -xM 
<m> 2 



(B.10) 



for all to G Z. The consequence of all this is that for indices ij € {0, 1, 2} and 
]3, 5=1, ...,77 — 1 one has 



I % (filial • • • \ik-l\VpVq) (TO) I < 7ii,i 2 ,-, 



g-x|m| 
<777> 2 



Vm G Z, 



from some positive constants 7.^. 



_ 1 , depending on the indices ij. 
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Returning our attention to expression (A. 18), we conclude that 



n 



(i) 

n+3 



< 7/2 2J K p K n ^ p + 7/3 ^ K p K n+1 - p + f]4 ^ KpKn+2- P , (B.ll) 

p— 1 p— 2 p— 3 

for certain positive constants 772, 773, 774. 

Bound for T 3 . Si nce qv 2 is a linear combination of the functions Q , Qx and 
Q2, we conclude from (A. 19) and from the previous arguments that 



n 



(2) 
n+3 



n-1 



< 775 KpK n -p + 116*^2 K p K ^+i-p ' 

p=l p=2 



for certain positive constants 775, ?yg. 

Bound for T 4 . By fl4.28|), one has 



(B.12) 



M(q 1 v 3 l n ) = i }^M(9j 2 (qv3 \v P v n - P ))- 



0=1 



From (4.21) and (B.10), we see that both qv% and v p v n - p satisfy the conditions 
of Proposition [si!! . Hence, by ( |B.S| ), M(9l2 (<Z W 3 \ v p v n-p)) < f]7KpK n - p for some 
positive constant 777 and 



M(q- X v z l n ) < r] 7 J2 K P K n-p- 
P =l 

We arc finished with the bounds for the terms T, of 



(B.13) 



If we collect flB.q), 

( B.ll ), ( B.12 ) and ( B.12 ) and return to (B^), we conclude that there are positive 
constants Ti, I^, I^, T4 such that we can recursively define 



K n := r 1 J2K p K n _ p + T 2 J2K p K n+1 _ 



P =i 



p=2 



+ r 3 2J K p K n+2 -p + r 4 2J KpK n+ 3- p , 
p— 3 p— 4 

for n > 4, after fixing the convenient values for i^i, K2, -K3 and K4. Note that 
we can choose K\ = K 2 = K3 = K4 taking Ki = max{K\, K 2 , K3, K4} for all 
i = 1, . . . , 4. Defining T := max{ri, T2, 1^, 1^}, we can redefine the K n so as 
to have the more convenient choice 



K n 



7i—l n—1 

KpK n _ p + k p k 

n+1— p 

_p=l p=2 

n—1 n—1 

+ ^2 KpK n+2 -p + 

p—3 p— 4 



(B.14) 
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Expression (B.14) has an analogous one in g, with the distinction that only 
the two first sums above occurred in the corresponding expression for K n . From 
now on, we follow closely the steps of 0. The first one is to show that K n is a 
non- decreasing sequence. We have 



K„ 



n+l 



L n+4— p ^n—p) 



n-l 

+ K p {K n 

p=i 

Therefore, assuming inductively K\ = K2 = K3 = K4 < . . . < K n implies K n < 
K n+ \, thus proving that the sequence is non-decreasing. Next, we write (B.14) as 



K n = T 



n-l 



12\12 Kb ) Kn - a + 12 k p ( Kn -p + Kn + 

i-p + K n+2 -p + K n+3 - p ) 



a=l \b=a 



p=4 



Since the sequence is non-decreasing, we have K n _ p +K n+ i- p +K n+ 2- p +K n+ 3_ p < 
4K n+ 3^ p and K n - a < -Kn-i for a = 1, 2, 3. Hence, we may say that 

3/3 \ n-l 
X « ^ r 12 [12 Kb ) Kn - 1 + 4r 12 K P K n+3- P 



a=l \b=a 



p=4 



+3- 



p ; 



(B.15) 



p=4 



where f := j Kf, j is a positive constant. Adding up the positive quan- 



a=l \b= 



71-2 



tity f K p K n+ z- p to (B.15) and setting A := max{f , 4r}, we get 

p=4 



K n < Aj2K p K n+3 _ p . (B.16) 

p=4 

Let us now define another auxiliary sequence Jk such that Ji = Ki for I = 
1, 2, 3, 4 and 

n-l 

J n . A ^ ^ JpJn+3—p • 
p=4 

for n > 4. It is a simple exercise to show from (B.16) that K n < J n for all n. Now, 
let us consider the translated sequence L n = J n +2, n > 1. We have 

n+l n+l n— 1 

L n = A ^""^ JpJn+5-p = A ^""^ L p -2L n +3-p = A ^""^ LpL n +i- p . (B.17) 



p— 4 



p=4 



p=2 
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n-1 

The sequence c„ defined by c„ = c p c„_ p+ i for n > 3, with Ci = C2 = 1, 

p=2 

defines the so-called "Catalan numbers" , which can be expressed in a closed form 

as 

(2n-4)! 

- - n > 2 . 



(n-l)!(n-2)!' 

By Stirling's formula, the c n 's have the following asymptotic behaviour: c„ rs 
16 ^P ^72 > f° r w large. The existence of a connection between the Catalan numbers 
and the sequence L n is evident from (B.17). Two distinctions are the factor A 
appearing in (B.17) and the fact that L\ = L 2 = K3 = K4 are not necessarily 
equal to 1. One can, however, easily show that L n = (i^)" -1 A"~ 2 c„, n > 2. 
Hence, the following asymptotic behaviour can be established: 

L ~ 1 ( 4 ^A)" 

Ln ~ 160F^ 3 A 2 n3/2 ' nlarge ' 

Since if„ < J„ = i n _2, we conclude that for n large K n < M (M 1 ) n , for 
some positive constants Mo, Mi. This completes the proof of Theorem ■ 



C Comments on the Fourier Coefficients 

Relations and can be obtained in our case by repeating the analysis 

of B. One of the conclusions is that the Fourier coefficients of ( |l.9| ) are analytic 
functions of e, for |e| small enough. There is, however, a point to be noticed here. 
If we follow the steps of |2j and compute the Fourier expansion of R(t)~ 2 to obtain 
the Fourier expansion for S(t), this last function will contain terms like J* e 2lQ,t dt, 
which behave like f2 _1 . In case (I), treated in j|, this is not problema tic, since 
O = O(e) and since such terms are always multiplied by a factor e (see ( |2.2| )). In 
our case, however, Q = 0(e 3 ) and we have to look such terms more carefully. 
Looking at the e-expansion for U\2(t), we have 

U 12 (t) = -ie~ 1 ^ [eWx(t) + e 2 W 2 (t) + 0(e 3 )] , 

where, after a lengthy computation, we get 

' i(nu) + 2S2) L 

n,m y ' 

The problematic term is that proportional to e/f2, which appears for n = 0. How- 

(2) 

ever, this term is proportional to Q — M(Qo) — 0, by hypothesis and, hence, 
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the Fourier coefficients of eW\ are analytic on e. Let us look now at e 2 W2- After 
another lengthy computation, one gets 

} o {2) o~o m 1 

_ e i(pu+n)t _ Wn-mW-pW™ i(pu-Sl)t 

m(nio + 20) 



(pw-O)t _ On Qm+n-pQm t(ph)+n)f 

m(nio + 2f2) 

Again, the dangerous terms are those with n = 0, since the are proportional to 
e 2 /f2. However, one sees by inspection that such terms vanish in the expression 

(2) 

above. In fact, in the last two sums they appear proportional to Qq , which is 

zero. In the first two sums, they appear proportional to X^m^o — m ~ ' which * s 
also zero, as one sees by interchanging m — > — m. 

Our conclusion is that all dangerous terms like I/O, appear multiplied at least 
by factors e 3 and are, therefore, analytic. 



|2 



2 W 2 (t) = 



E 



Qn—mQn—pQ- 

m(nuj + 2f2) 



— T. 

p, n, m 



luJ 



o {2) o — o {2) 

m(nu> + 2f2) 



D An Identity on Sums of Bessel Functions 



Here we will sketch the proof of identity (5.2), since we did not found mention to 
it in the literature. Using the well know identity (due to Schlafli and Gegenbauer) 
for products of Bessel functions (see frflj ) 



i r' 2 

J ll {z)J v {z) = - Jfj, +l/ (2zcos9)cos((n- v)6)d0, 

T J -tt/2 



(D.l) 



one gets 



S m (x) := 

k^O 



J m+k (^-) 



(-1)' 



tt/2 



-tt/2 



Jo{2x cos 9)g(9)e 2mU d$, 



where g(9) := ^^(-l)'^. On the interval (-tt/2, tt/2) one has g{9) = 
-26 /w. Hence, 



2(-l) m 






-7t/2 




(a 


7T 





J (2a; cos 6»)6> sin(2m6>) d6» 

7T/2 



-tt/2 



J (2x cos6») cos(2i/6») d6> 



Using again ( D.l ) (with /x = —u), one finds 



S m (x) = (-1) 



[J v {x)J- v {x)} 



J in (#) 



-2^-J v (x) 
av 



+ nY m (x) 
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This is Q5.2D, as desired, where Y m is the Bessel function of second kind, given by 
= UkM*) ~ (- l ) m l J ^))l^ for integer m. 
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